Some remarks on Oscillating Inflation 



V. Cardenas* and G. Palma^ 
Departamento de Fisica, Universidad de Santiago 
Casilla 307, Correo 2, Santiago, Chile. 
(February 1, 2008) 

In a recent paper Damour and Mukhanov describe a scenario where inflation may continue during 
the oscillatory phase. This effect is possible because the scalar field spends a significant fraction 
of each period of oscillation on the upper part of the potential. Such additional period of inflation 
could push perturbations after the slow roll regime to observable scales. Although in this work we 
show that the small region of the Damour-Mukhanov parameter q gives the main contribution to 
oscillating inflation, it was not satisfactory understood until now. Furthermore, it gives an expression 
for the energy density spectrum of perturbations, which is well behaved in the whole physical range 
of q . 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Nowadays inflation is a widely accepted element of 
the early cosmology . It gives the possibility of solv- 
ing many of the shortcomings of the standard hot big 
bang model and provides the source for the early energy 
density fluctuations responsible of the large scale struc- 
ture of the universe observed today. Although there are 
many models of inflation, the underlying physical ideas 
are well established. These are characterized by a period 
of "slow roll" evolution of a scalar field (called inflaton) 
toward the vacuum potential. During this period the 
field changes very slowly, so that the kinetic energy ip 2 /2 
remains smaller than its potential energy V(ip). The en- 
ergy density associated to the scalar field acts as a "cos- 
mological constant" term, allowing a period of quasi ex- 
ponential expansion of the scale factor. When the period 
of inflation ends, the scalar field tp start a phase of rapid 
coherent oscillations around the vacuum. 

Very recently ( ]^||) it nas been pointed out that in- 
flation can persist during the coherent oscillations of the 
inflation field phase. This exciting result is possible when 
the inflaton potential verifies a simple constrain of curva- 
ture far from the core convex part, where the inflaton field 
can roll slowly. The efficiency of this phenomena could 
have important implications for GUT scale baryogene- 
sis ill. In fact, as suggested by Damour and Mukhanov 
( j,it can be expected that due to the increase of the 
oscillation frequency, there is the possibility to generate 
massive particles heavier than ~ 10 16 GeV. 

In ref. Q Damour and Mukhanov estimated the 
amount of inflation to be ~ 10 e-fold (powers of the scale 
factor). They argue that this effect can be more efficient 
than the parametric resonance effect |a] for the amplifi- 
cation of cosmological perturbations |6 1 . In ref. || Liddle 
and Mazumdar showed that Mukhanov et al. overesti- 
mated the number of e-fold because they have used a 
slow-roll definition of this object. In their paper, Liddle 
and Mazumdar found an analytical expression for the 



number of e-fold of inflation using the appropriate defi- 
nition finding a number of ~ 3 e-fold concluding that this 
effect is not very efficient. The study of adiabatic pertur- 
bations in this phase has been made by Taruya [Q. He 
found a poor amplification in the case of a single scalar 
field model but anticipated an enormous amplification for 
multi-field systems. 

In this letter we review the problem. In particular 
we find that the analytical expressions used to compare 
with the numerical estimation are not well defined in the 
q ~ region and propose a way to correct these analytical 
estimations. Furthermore, with this result we study the 
evolution of the scalar field finding total agreement with 
the conclusions of ref. || for q > 0.2, but a remarkable 
different result for small q. For this region, the initial 
conditions are very important. We find that q ~ gives 
the leading contribution for oscillating inflation and the 
dominant part in the amplification of the fluctuations. 

The letter is organized as follow; first we describe 
briefly the Damour-Mukhanov model. Then, we make 
some comments about the initial conditions for this phe- 
nomenon and later we propose an improved expression, 
valid for the leading region of q, which is our main con- 
tribution. 



II. BASIC EQUATIONS 

Now we shall restrict ourselves to models of inflation 
driven by a single scalar field. The equations are 



(p + mtp + v ilp = o, 



(1) 



(2) 



Here H = a/ a is the Hubble parameter, a is the scale 
factor of the universe and k 2 = 8ir/3M 2 with M p — 
1.2 • 10 19 Gev the Planck mass. During the oscillatory 
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phase of <p we have two time scales; the inverse of the 
frequency w _1 of oscillations of ip and the inverse of the 
rate of expansion H^ 1 . If the limit ui ^> H is taken we 
can neglect terms proportional to H in the equations. So 
from (1) we can integrate to obtain, 



1 -2 



V = cte = Vrr, 



(3) 



where V m — V{<p m ) is the maximum value of V(ip) in 
each oscillation when the field reaches the maximum 
value tp m . From this relation we obtain the period of 
a single oscillation, T = 4 f* m dip [2(V m - V(<p)]* .When 
lu <C H we can define an adiabatic average index 7 by 7 = 
(G° + P)/p) , where the bracket means (...) = T -1 ...dt. 
Equations (1,2) can be re-written in the fluid form 



-3H(p + p), 



a 1 
a 3 



•3p), 



(4) 



(5) 



then from the definition of 7 and eqns.(^,|^) we have sev- 
eral ways to compute the adiabatic index 
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<^ 2 ) (<pV, v ) 



Vrr, 



Vrr, 



(6) 



Because p = (7 — l)p and (g) we have a super luminal ex- 
pansion a > when 7 < 2/3. From the last two relations 
in eqn.(|6|) the inequality 7 < 2/3 leads to 



{V-ipV, v ) >o. 



(7) 



III. THE DAMOUR-MUKHANOV MODEL 

Until now everything has been done for an arbitrary 
potential, but from now on we shall consider the potential 



V(<p) 



A 

'1 



9/2 



(8) 



where q is a dimensionless parameter, A = [mass] 4 is 
a constant and ip c = [mass] determines the size of the 
convex core of V(ip). We assume for a while that tp c 
marks the end of oscillating inflation. The analysis made 
in ref. Q works well far from the core of the potential. 
Further, the limit <p 3> <p c of eqn.(H) was written as: 



V(<p) 



A ftp}* 
q \<Pc 



(9) 



In this case, the adiabatic index can be computed exactly 
given || by 
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_2q_ 

g + 2' 



(10) 



so, from the inequality 7 < 2/3 we note that to hold 
inflation during the oscillatory phase we must have q < 1. 
By using eqn.(jiTj|) in eqn.(Q) we obtain p = —3H~fp and 
together with eqn.(|J) we have 



a oc t 2 '^ = ^+ 2 )/ 3 ? 



ipm OC t 



■2/9 « -6/(fff2) 



(11) 



(12) 



(13) 



where ip m is the amplitude of the oscillations, ip c < tp m < 
ip s and ip s is a typical value of co at the end of slow-roll 
inflation and the beginning of oscillating inflation. To 
compute the number of e-fold of inflation during oscil- 
lating inflation we can not use the standard expression 
N = ln(a//ai), appropriate for the slow-roll stage, but 
the improved expression proposed in ref. H 



N = In 



a f H f 
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(14) 



because in each oscillation, while the field spends time in 
the core region, the universe continue their expansion so 
H can vary. Then from (§) and @ H oc a -3 9 /( 9 +2) the 

product aH oc tp, 



(i-g)/3 



and from f rL4|) we obtain 



N 



In 



qM p 



(15) 



where we have used ip s ~ qMp/yWir. In |3| the numer- 
ical curves for co c = 10~ 6 -M„ show that iV < 3. Using 
the analytical expression ([l5|) we do not find agreement 
for small values of q. However there is not a compelling 
reason to believe in (Il5|) for small q. 



IV. THE SMALL Q-REGION 

To study the small q region, we must use the correct 
limit q — > of (||) which leads to 



£| +1 



(16) 



so, if now we take the limit p S> p c we obtain the loga- 
rithmic potential VUp) ~ Aln(cp/cp c ). A very important 
fact to note from (Q) is that the limit c/ — > does not 
exist. Of course, the expression (^J) is wrong around the 
q ~ region and the expressions derived from this are 
ill-defined. But, some work has been done in this regard 
For the logarithmic potential the adiabatic index is 
7 = l/\n((p m /(p c ), so from (Q) and (|2|) we obtain 



a(t) oc exp 



/I 



{tend — t) 2 



(17) 
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but this form does not permit us to write an explicit ex- 
pression for N (see ref . |3j] ) . Let us make some comments 
about this result. Because 7 = 1/ \\i{tp m /ip c ) from (||) 
we obtain p = —3HA, then a ~ (^m/^Pc)^ 1 ^ 3 ■ More- 
over, from (||) we have H oc p 1 / 2 oc (lnf^m/^c)) 1 / 2 , 
then to compute N we should evaluate the factor 
(ln(v? m /( ( c c )) 1 / 2 ((/3 m /( / 3 c ) _1 / 3 at the extremes ip s and tp c , 
but this is not possible. The problem arises when ip m 
is chosen close to tp c in a expression valid for tp ip c . 
Because ([)[) is valid for ip ^> ip c too, the same problem 
should appear in the calculation of N. In fact this is the 
case but, to see that, we must include the constant term 
A/q in (|9|) or equivalently, write the limit correctly. 
When we use the potential 



^\ -1 



(18) 



the adiabatic index become time-dependent, satisfying 
the equation 



2q A 



(19) 



where p(t) = V{w m (t)) = A[(tp m (t)/tp c y - l]/q. Re- 
placing this in (0) we obtain the same behavior as in 
eqn. (|l|). But when we calculate H we obtain H oc 
p 1 / 2 oc q~ 1 / 2 [{ip m /(p c ) q — l] 1 / 2 , which is not well defined 
at the point ip m = (p c . The same happens to the e-fold 
number, as mentioned before. 

The authors of ||] found ip s ~ qAl p /yfl6ir using the 
potential (^). From this result they found a decrease of 
N at decreasing values of q. This seems very strange 
because, for smaller values of q we obtain flatter poten- 
tials at (p 3> (p c , then tp rolls slowly most of the time 
increasing the amount of inflation. Furthermore, we are 
not safe of how they set the initial conditions for ip in 
their numerical analysis. In general, the initial and final 
field configuration depends on the potential. 

For q close to cero, tp s is not proportional to q. In fact, 
we know that tp s came from the saturation of the slow- 
roll inequality \V /V\ < \J4&-k/M p using the expression 
( |l8| ) for the potential 



q-l 



qM P 



(20) 



for the ip ^> ip c region. In Figure 1 we see the behavior 
of ip s in terms of q given by eqn.(po|). Of course, in the 
large g-region both curves agree. In order to illustrate 
this point and compare with ref. |^] we use its expression 
given by eqn.([lj) but instead of using ip s ~ qM p /\/T6ir 



we use the numerical values of ip> s obtained from eqn. (20). 
The results are plotted in Figure 2. In the small g-region 
the field <p s grows preventing the fall of N predicted in 
ref. §. 

Moreover, as we have anticipated before, the value of 
the field at the end of oscillating inflation will have a q- 
dependence too. We know from ref. j| that the intercept 



U(ip) = V(ip) — (pV. v , must be positive to hold oscillating 
inflation. Let us define <p>f to be the value of the inflaton 
field <p at which U(<pf) = 0. This condition represents 
the end of inflation due to oscillation. We need thus to 
compare ipj for different values of q. 

If we take the potential ( |l8| ) and define x = <p/tp c , we 
obtain 



U(x) = - [X« (1 - q) - 1] . 



(21) 



From this equation we can extract a explicit expression 
for ipf. If we impose U(xf) = we obtain the value of 
the scalar inflaton field at the end of this phase 



iff =tp c (l- q) 



-l/g 



(22) 



Using the improved expression eqn. (0) for the e-fold 
number (ref. |3|), but inserting the corrected values for 
tp s and ipf given by eqns. ( |20| ) and ( |22| ) we obtain 




(2+ 9 )/6 



- 1 

LW<^) 9 -i 



1/2^ 



(23) 



The corrected value for ip f leads to a even smaller amount 
of inflation when comparing with the value obtained in 
ref. J3|. In figure 3, we plot eqn.(|23|) and show the be- 
havior of both effects combined. Because ipf is greater 
than tp c , the amount of inflation is smaller than one ob- 
tained by Liddle et al. || in the whole range of q € (0, 1). 
Moreover, the correct values of <p s produce a positive con- 
tribution to N in the small g-range, which avoids the fall 
of N, as q goes to smaller values, predicted by Liddle et. 
al. ||. Again, it is not possible to show the whole range 
of q because q > (p c /v ( see comments below eqn.(po|)). 

Because oscillating inflation adds e-folds of inflation 
after the slow-roll regime, where the observed perturba- 
tions are generated, is possible that this additional period 
of inflation could push perturbations to observable scales. 
In order to obtain the required amplitude of density per- 
turbations and without imposing unphysical constraint 
on the potential, we should compute the density pertur- 
bation spectrum for the model being studied. 

In reference || an expression for this object was de- 
rived 



51 = 



512tt A 



9+2 



75 q 3 Mji tpl 



However, this expression is not well defined close enough 
to zero. Using the primordial density perturbation spec- 
trum Sh as was define in M we obtain for eqn. 



51 = 



512tt A 



75 



1 3M P 1 (<p/<Pc) 



2q 



which is well defined even for the small values of q 



512tt A 
75 Mjj 



ln'1^. 

•Pc 



(24) 



(25) 
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Because the COBE satellite require Sh ~ 2 • 10~ 5 , 
in the q — case the amplitude of the potential for 
ip c = l(T 6 Mp gives A 1 / 4 ~ 2 ■ 1CT 3 M P , which is a typical 
number for inflationary models. 

V. SUMMARY 

As a summary, we have made some corrections about 
how to compute the e-fold number, which accounts for 
the amount of inflation during the oscillatory phase. In 
particular we note that previous studies are not accu- 
rate because they are not valid close to the core of the 
potential ip ~ ip c . Thus, we make the analysis for the 
small q region of the potential, which has not been consid- 
ered until now. Finally we find that, in order to extract 
the correct amount of inflation during this phase, a very 
careful definition of initial and final field configuration 
is needed. Our results show that near q ~ the e-fold 
number is maximal but it is still not enough to be more 
efficient than the parametric resonant effect discussed in 

!■ 



VI. FIGURES CAPTION 

Figure 1: We plot the q dependence of ip s . We see 
that for q > 0.2 both curves agree but for smaller values 
of q the field ip s grows, preventing the fall of N. 

Figure 2: The e-fold number N is shown as a func- 
tion of q, taking into account the behavior of ip s described 
for the eqn. (20). 

Figure 3: The e-fold number N is plotted vs q, 
taking into account the combined effects: the behavior 
of ip s described by eqn. (20) and the definition of iff 
discussed in the text. 
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